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Abstract. The Kustaanheimo-Stiefel (KS) transformation maps the non-linear 
and singular equations of motion of the three-dimensional Kepler problem to the 
linear and regular equations of a four-dimensional harmonic oscillator. It is used 
extensively in studies of the perturbed Kepler problem in celestial mechanics 
and atomic physics. In contrast to the conventional matrix-based approach, 
the formulation of the KS transformation in the language of geometric Clifford 
algebra offers the advantages of a clearer geometrical interpretation and greater 
computational simplicity. It is demonstrated that the geometric algebra formalism 
can readily be used to derive a Lagrangian and Hamiltonian description of the 
KS dynamics in arbitrary static electromagnetic fields. For orbits starting at the 
Coulomb centre, initial conditions are derived and a framework is set up that 
allows a discussion of the stability of these orbits. 



PACS numbers: 45.20.Jj,31.25.Gy,45.50.Pk 

1. Introduction 

The Kepler problem belongs to the simplest systems of classical mechanics. A detailed 
understanding of its properties is of equally fundamental importance to celestial 
mechanics and to atomic physics. In both areas it is of interest to describe, either by 
means of analytic approximations or numerical computation, the impact of additional 
non-Coulombic forces on the dynamics. In its original form, the equation of motion 

d 2 x x 

d^ = ~]xj ( > 

is not well suited to this purpose because it is highly non-linear and exhibits a 
singularity at the Coulomb centre where the force diverges. For numerical studies 
of the dynamics it is mandatory to find a representation of the equations of motion 
which avoids this singularity. 

For the one-dimensional Kepler motion, it was already found by Euler 1 that 
the introduction of a square-root coordinate u — \/x and a fictitious time r defined 
by dt — x dr reduces the Kepler equation of motion Q to the equation of motion of 
a one-dimensional harmonic oscillator 
d 2 u 

— +2Eu = 0, (2) 
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where E is the energy of the Kepler motion. Equation @ is not only void of 
singularities, it is also linear and thus forms a much more convenient basis for analytic 
calculations. 

Generalizing this approach, Levi-Civita [2] regularized the two-dimensional 
Kepler problem by combining the two spatial coordinates into a complex number 
x = X1+1X2 and introducing a complex square-root coordinate u = y/x, which together 
with the fictitious-time transformation dt = \x\ dr reduces the Kepler problem to a 
two-dimensional harmonic oscillator. 

Attempts to extend this regularization scheme to the three-dimensional Kepler 
problem failed, until in 1964 Kustaanheimo and Stiefel |3I4| proposed the introduction 
of four regularizing coordinates instead of three and thereby achieved the reduction of 
the three-dimensional Kepler problem to a four-dimensional harmonic oscillator. This 
transformation, which is known as the Kustaanheimo-Stiefel (KS) transformation, is 
discussed in detail in the monograph by Stiefel and Scheifele : 5 . Beyond its importance 
to celectial mechanics, it has proven to be an essential tool for investigating the 
complicated classical dynamics of the hydrogen atom in crossed electric and magnetic 
fields 0C1IH1- 

Customarily, the KS transformation is expressed in the language of matrix 
algebra, which not only necessitates awkward computations of vector and matrix 
components, but also lacks a transparent geometric interpretation. An alternative 
formulation in terms of the geometric algebra of Euclidean three-space was introduced 
by Hestenes In this formalism, the four KS coordinates are interpreted as the 
components of a position spinor and are thus given a clear geometric meaning. In 
addition, the formalism offers computational advantages over the conventional matrix- 
based approach because it unites the four coordinates into a single spinor. 

In this paper I elaborate on the geometric algebra formulation of the KS 
transformation and work out the details necessary for applications to atomic physics, 
in particular semiclassical closed-orbit theory [l()lllllH21ll3H14lll5| . The calculations 
will amply demonstrate the advantages of the geometric algebra over the matrix- 
based formalism. The terminology used is adapted to applications in atomic physics. 
In particular, I refer to the attractive centre as the nucleus and discuss the motion of 
an electron of unit mass and negative unit charge under the influence of the nucleus 
and arbitrary static external electromagnetic fields. Nevertheless, the results are of 
much broader range. They apply equally to celestial mechanics or any other field of 
physics governed by the equations of motion of a perturbed Kepler problem. 

In section [21 the geometric algebra formulation of the KS transformation is 
introduced and its relation to the matrix-based approch is established. The spinor 
equation of motion given by Hestenes [H] is derived. In section a Lagrangian 
and Hamiltonian formulation of the KS spinor dynamics is derived. It cannot be 
obtained by a simple change of variables because the KS transformation introduces 
a non-physical fourth degree of freedom and a pseudotime parameter. Nevertheless, 
it is demonstrated that the geometric algebra formulation readily lends itself to an 
incorporation into the Lagrangian and Hamiltonian formulations of dynamics. At 
the same time, the well-known KS Hamiltonian |12lllHl l6*]. which is restricted to 
homogeneous external electromagnetic fields, is generalized to arbitrary static fields. 
Section 0] presents the explicit solution of the spinor equation of motion in a pure 
Coulomb field and gives the constants of motion in terms of the KS spinor. Section [S] 
discusses the KS description of orbits starting at the Coulomb centre and returning to 
it. These orbits, referred to as "closed orbits", are of particular importance because 
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they play a central role in the semiclassical interpretation of atomic photo-absorption 
spectra pT71ITTlfT51fT^| . They require a special treatment because, although the spinor 
equation of motion is regular at the nucleus, the KS transform is singular there. I 
will derive initial conditions for orbits starting at the nucleus and present a basis of 
the spinor space that is suitable for the investigation of the stability of closed orbits. 
A brief introduction into the properties of geometric algebra needed here is given 
in |Appendix A| where the notation used in what follows is also explained. A more 
detailed exposition of the formalism can be found in [??l ll()l[T7| ll8j . 



2. The spinor equation of motion 



The KS-transformation in three dimensions can be found by representing an arbitrary 
position in space not by its position vector x, but by a position spinor, i.e. the 
rotation-dilatation operator transforming a fixed reference vector into the position 
vector x. As explained in |Appendix A| a rotation-dilatation of the reference vector 
<t 3 is represented in the geometric algebra by an even multivector U according to 

x = ~ Ua 3 U^ . (3) 

The factor 1/2 was introduced here to stay in touch with earlier applications of the 
KS-transformation to atomic dynamics [121113) . although the present formulation of 
the theory would suggest dropping it. It implies the normalization 

U^U = C/C/ 1 = 2r = 2\x\ . (4) 

Up to normalization, the ansatz J3J) reproduces the square-root coordinates introduced 
by Euler and Levi-Civita, respectively, if it is applied to spaces of one or two 
dimensions. 

Given a position vector x, the choice of the spinor U is not unique. More precisely, 
the gauge transformation 

U i — ► Uc- Isa/2 (5) 

with arbitrary real a does not alter x, because the additional exponential factor 
describes a rotation of the reference vector cr 3 around itself. This consideration 
immediately clarifies why a position spinor representation in three dimensions must 
introduce a fourth degree of freedom. In lower dimensions, a rotation does not leave 
any vector invariant, so that the spinor transformation does not possess a gauge degree 
of freedom. It is also clear from J5J that all fibres of the KS transformation except 
[7 = are circles in spinor space. 

The inverse KS transformation can be found from equation (|A. 13)1 by adapting 
the normalization to The position spinors corresponding to a vector x are given 

by 

u= r + x*s e _ l3a/2 
\Jr + z 

with arbitrary real a. 

In components, the spinor U can be represented as U — uq + Iu with u = 
X)fc=i u k (T k- The transformation then decomposes into 



X — UlU 3 — U U 2 , 

y = uiu + u 2 u 3 , 

1 / 2 2 2 i 2\ 
z = TjK - u l ~«2 + u 3> 



(7) 
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Up to renumbering the components, this agrees with the conventions of [121113) . 

To obtain an equation of motion for U, time derivatives of U must be calculated. 
Differentiating © leads to 

x = ~(r*3U 1 + \u<r 3 tf = {Utr 3 U i )i. (8) 

Equation J3J obviously cannot be solved for U because the time derivative of the gauge 
parameter a in JSJ cannot be determined from the dynamics of the position vector. 
To arrive at an equation of motion for U, I must therefore impose a constraint on a. 
This can be done in a convenient and geometrically appealing fashion by requiring 

{U<r 3 U% = , (9) 

which means that U is chosen such as not to contain a component of rotation around 
the instantaneous position vector x. Under this constraint, (jSJ) yields 

x = Ua 3 tf (10) 

and 

U = xU^ 1 <t 3 =x^-<t 3 . (11) 
2r 

As in the one- and two-dimensional cases, the regularization of the three- 
dimensional Kepler motion requires the introduction of a fictitious-time parameter 
t. It is defined by 

dt = 2rdr (12) 

Derivatives with respect to r will be denoted with a prime. Equation (|ll|l then yields 

U' = 2rU = xUa 3 . (13) 
For the second derivative of U, I obtain 

U" = (4p<k) Uct 3 + xU'ct 3 

= 2rx U(T 3 ^- U + x 2 U (14) 
2r 



/.. 1 



= 2 \^xx + -x'j U . 

Together with Newton's equation of motion 

x = ~ + f (15) 

with an arbitrary non-Coulombic force /, equation (|14|l yields the spinor equation of 
motion in the form first given by Hestenes 9 : 

U" = 2(E K + fx)U , (16) 

where the Kepler energy 

Ek = \* 2 -- (17) 
2 r 

denotes the sum of the kinetic and Coulombic potential energies. 

In the special case of pure Kepler motion, i.e. / = 0, the Kepler energy Ek is 
equal to the total energy E and is conserved. In this case, 1)16(1 reduces to the linear 
equation of motion 

U" = 2EU . (18) 
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If E < 0, this is the equation of motion of a four-dimensional isotropic harmonic 
oscillator with frequency u> = \ / —2E with respect to t. 

If additional forces / are present, the Kepler energy is not conserved in general, 
so that the work done by the external forces must be taken into account jjj. This 
can easily be achieved if the external forces are generated by static electromagnetic 
fields, because the work done by a magnetic field B is zero, whereas an electric field 
F = — W can be derived from a potential V(x). In this case, the total energy 
E = Ek — V is conserved, so that the equation of motion reads 

U" = 2(E + V{x) + fx)U (19) 

with / = —F — x x B. 

It finally remains to verify that the equation of motion is consistent with the 
constraint @. To prove this, first note that 

£ = {U'*zU% (20) 

is a constant of motion for any external forces / |3], because by Hlfijl 

^ = (U"* 3 tf) 3 + (U'a 3 U% 
=0 by (|XT2)l 

= ((E K + f x)x) 3 ( 21 ) 
= (E K x + r 2 f) 3 
= . 

Therefore, if the initial conditions are chosen so that £ = at r = 0, equation ( 1211 
guarantees (U cr 3 U^) 3 = 2r£; = at all times. 



3. Canonical formalism 



In classical investigations of atoms in external fields, the Hamiltonian nature of the 
dynamics plays a central role. It is therefore essential to show how the spinor equation 
of motion found in the previous section can be derived in the context of a Lagrangian or 
Hamiltonian formalism. In the matrix theory of the KS transformation, a Hamiltonian 
formulation is well known and widely applied in the literature |12II13| . Due to the 
introduction of an additional degree of freedom and a fictitious-time parameter, it 
cannot be found by a straightforward change of variables. In this section it will be 
shown that the application of geometric algebra allows an easy and general derivation 
of the Hamiltonian. At the same time, the Hamiltonian formalism will be generalized 
to arbitrary inhomogeneous static external fields. 



3.1. Fictitious-time transformations 

Elementary expositions of Lagrangian and Hamiltonian dynamics usually treat the 
time t as the externally prescribed independent variable fundamentally different from 
the spatial coordinates, velocities, and momenta. The formalisms are then shown to 
be invariant under point transformations or canonical transformations, respectively, 
which may be time-dependent, but may not transform the time variable. However, 
both the Lagrangian and the Hamiltonian formalisms can be reformulated in such 
a way that it is possible to introduce an arbitrary orbital parameter r and to treat 
the physical time t as an additional coordinate on the same footing as the spatial 
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coordinates. This formalism is discussed in its full generality by Dirac |19U2U| . For 
the special case of autonomous Lagrangian dynamics and the simple form of the 
fictitious-time transformation used above, the full flexibility of Dirac's homogeneous 
formalism is not needed. Instead, the modifications needed to achieve the fictitious- 
time transformation can be derived in a straightforward manner from Hamilton's 
variational principle. 

The Lagrangian equations of motion can be derived from the action functional 

S = f " dtL(q(t),q(t)) (22) 

by requiring that for the classical paths the variation of S with respect to the path q(t) 
vanishes if the variation is performed with the initial and final times t% and £2 and the 
coordinates q(t\) and qfa) kept fixed. If a fictitious-time parameter r is introduced 
by the prescription 

dt = f{q,q)dr (23) 

with an arbitrary function /, it is tempting to rewrite the action functional as 

S = f 2 drf(q,q)L(q, q ) (24) 



and regard 



L = f(q,q)L(q,q) = ^-L (25) 

as the Lagrangian describing the dynamics with respect to r. However, this simple 
procedure is incorrect in general, because to derive the Lagrangian equations with 
respect to L from l|24|l . the variation of S has to be performed with the initial and 
final fictitious times T\ and T2 kept fixed, and due to (|23|l a variation of the path will 
alter the relation between t and r, so that the initial and final physical times t\ and 
t2 will vary. 

To establish the true relation between (12211 and l|24JI . I calculate the variation of 
H24[l taking the variation of t into account, i.e. q and t are varied according to 

q(r) i-> q(r) + 5q(r) , 

i(r) 1 — ^ *(t) + 6t{r) , 
subject to the boundary conditions 

5q( n ) = 6q(T 2 ) = (27) 
and with n and t 2 kept fixed. Under this variation, 

dr clr 1 

dt ^ d(t + St) = M + p 

v ' dr dT (28) 
dr ( dr di5t N 



dt \ dt dr 



so that 



dr / dr \ 2 d St 



dt V dt J dr y ' 
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and hence 



dq \ dr dq 

,dx/ ' d7 + dr" 
d(5q . dr d 5t 

= ~dT 

The variation of i(24|) then reads 



(30) 



q dt dr 



55 = / dr 



L+ — 



dt 

d7 



dr — (l - ' — 
dr 1 ^ <9a 



— 5g- 
dq 



dt 



_dL s . 
dq 

dL § + 8LdSq 
dq dq dt 



(31) 



dt 



dL _d_ dL s s 
dq dt dq 



where the customary partial integration was performed, the boundary conditions (|27|l 
were used and the Hamiltonian 

H = q^-L (32) 

was introduced. 

If only the second integral in the last line of 1(31(1 were present, it would yield the 
correct equations of motion. Thus, the action functionals 1(22(1 and 1(24(1 are equivalent 
if the Hamiltonian H vanishes. For autonomous systems, H is a constant of motion 
equal to the energy E. If the Lagrangian L is replaced with L+E, with E regarded as a 
constant, the equations of motion derived from L are unchanged, but the Hamiltonian 
1(32(1 changes to H — E = 0. Thus, the dynamics of trajectories with energy E with 
respect to the fictitious-time parameter r is described by the Lagrangian 



C=^{L + E) = f{q,q) (L + E) 



(33) 



This Lagrangian has to be written as a function of the coordinates q and the fictitious- 
time velocities q' . If the function / is independent of the velocities, the canonical 
momenta are invariant under the fictitious-time transformation, because q' = f(q) q 
and 



dC _ dq dL _ dL 
dq 1 ' dq' dq dq 



(34) 



From the time-transformed Lagrangian 1)33(1 , the transformed Hamiltonian 

f(q,q)(H-E) (35) 



H 



is obtained by the usual Legendre transformation. It must be written as a function 
of the coordinates and momenta. In some cases the passage from the Lagrangian 
to the Hamiltonian description of the dynamics is impossible because the relation 
p = dC(q,q')/dq' cannot be solved for q 1 . In these cases, the Hamiltonian 1(35(1 can 
be shown to describe the fictitious-time dynamics by a discussion of the modified 
Hamilton's principle analogous to the derivation of the Lagrangian C above. 
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3.2. Lagrangian description 

The dynamics of an atomic electron under the combined influences of the nuclear 
Coulomb potential, an additional scalar potential V(x) and a magnetic field 
represented by a vector potential A(x) is described by the Lagrangian 

x 2 1 

L = — + - + V(x) - A x . (36) 
2 r 

This Lagrangian must be transformed to a Lagrangian £ describing the fictitious-time 
dynamics of the position spinor U. With f(q) = 2r = WU by i|12|) . the fictitious-time 
Lagrangian (|33|l reads 



C 4 



- U*U' + — / (U'<t 3 U^) 2 ) + EU^U + U^U V(x) - (A(x) U'a a rf) + 2 (37) 
4 8r \ I 

with x — ^Ucr 3 U' f . If the constraint © is used, £ simplifies to 

£ = - ET'V + £ J7 f J7 + C/ f C/ V(x) - (A(x) U'(t z U ] ) + 2 (38) 

Both forms of the Lagrangian yield the same "on-shell" dynamics for trajectories 
satisfying ©. Note that only the kinetic term is influenced by the constraint, whereas 
potential and vector potential terms are not. 
The momentum conjugate to U is given by 

P = du>£ = ^- trarfiU'asrth - ^A , (39) 
2r 

which simplifies to 

P = U' f - (T 3 U^A (40) 

if J3J is applied. 

As the spinor equation of motion l|16() is valid under the constraint © only, the 
Lagrangian £ provides a suitable description of the dynamics if it reproduces (jltjfl 
for trajectories satisfying 10. The simplified Lagrangian (|38|l can therefore be used. 
When equations of motion are derived from l|38|) , the constraint 10 must be taken into 
account by a Lagrangian multiplier. I will now show, however, that the unconstrained 
equation of motion 17 

^-d u ,£-d u £ = (41) 

derived from l|38|l reproduces Ijltjfl without the constraint being explicitly dealt with, 
i.e. the Lagrangian multiplier to be introduced turns out to vanish identically. I 
therefore ignore it from the outset. 

For the case of vanishing external potentials, (|41|l can easily be seen to yield 

U" j - 2EU^ = , (42) 

which is the reversion of H18|) . For the terms containing the potentials, the calculation 
of H41|) is still straightforward, but requires a more intimate familiarity with the 
properties of the multivector derivative. I will therefore present the calculation in 
detail. 

The contribution of the scalar potential term 

V = U ] UV{x{U)) (43) 
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with x(U) = \Ua 3 W reads, by ipO^I . 

d v V = 2U ] V{x) + tfUduV{x{U)) . (44) 
The chain rule l)A.24(l then yields for any even multivector M 

M * duV(x(U)) = (M * du x(U)) * d x V (45) 

with 

M * d v x(U) = (M * du)- Uct 3 U ] 
= (Ma a rf)i . 

In the absence of a magnetic field the external force is / = d x V, so that 

M * duV(x(U)) = ((MtraC/t),/) = (M* 3 tf f) . (47) 

Thus, 

duV = d M (M * duV) = crstf f , (48) 

and finally 

duV = 2U i V(x) + 2U*xf . (49) 

This is the reversion of the scalar-potential terms in Therefore, l|49ll together 

with (|42() indeed yields the correct equation of motion. 

To evaluate the contribution of the vector potential term 

A= (A(x)U'(t 3 U^ (50) 

to 14111 . first note that 

WW) (5i) 

= <r 3 U* A(x) + a 3 U\x' ■ d x )A(x) . 

By Leibniz' rule and ©, 

duA =du (A U'cr 3 du (A U'<r s rf) (g2) 

= a 3 U' j A + du{A(x(U))x') . 

The second term on the right-hand side can be evaluated by first calculating directional 
derivatives. For an arbitrary even M, (|46|) and the chain rule i|A.24|) yield 

(M * du)(A(x(U)) x 1 ) = (M * d v A) * d A (Ax') 

= ((M *duA{x(U)))x') 

= {{M * d v x) * d x Ax') (53) 

= (((MersU^^Ax') 

= (M(T 3 U^d x )(Ax') , 

so that 

du(A(x(U))x') = d M (M*du)(A(x(U))x') 

* , ( 54 ) 

= a 3 U^ d x (A ■ x') . 
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Equations J3TJ , and combine to 

4- d v ,A - duA = a 3 tf ((as' • - fiU(A • *')) 

= -<T 3 [/ f x (a. x A)) ( 55 ) 
= -2C/ f a;(i; x B) , 

which is the reversion of the magnetic-field contribution to 1)160. 

Note that the derivation given here is valid for arbitrary external potentials 
V and A, whereas conventional treatments restrict themselves to the special case 
of homogeneous external fields. It can obviously be further generalized to include 
conservative forces of non-electromagnetic origin. Also note that the geometric algebra 
formalism allows one to do the calculations in a straightforward manner without having 
to resort to component decompositions of any of the vectorial or spinorial quantities 
involved. 



3.3. Hamiltonian description 

The transition from a Lagrangian to a Hamiltonian description of the dynamics leads 
from the Lagrangian i|37|) or Q3H|). depending on whether or not the constraint © is 
applied, to the Hamiltonian 

H = (U'*dui)£-£, (56) 

in which the velocity V has to be expressed in terms of the momentum P. The 
transformation requires that the relation l|39l) or (|40|l between velocity and momentum 
can be solved for the velocity, which is impossible in the case of (|39J) . Thus, the 
constraint J^J) is not only needed to obtain an unambiguous equation of motion for U, 
but also serves as a condition for a Hamiltonian description of the spinor dynamics to 
exist. If it is imposed and an inessential constant of 2 is added, the Hamiltonian reads 

H = - (P f + A Ua 3 ) (P + o- 3 [/ f A) - EU*U - U^U V{x) = 2 . (57) 

Because it is time-independent, the Hamiltonian (|57|l is a constant of motion. To 
describe the physical dynamics, its value must be chosen to be 2, whereas the physical 
energy E appears as a parameter in TL. 

The equations of motion derived from l|57(l read 

U' = d P H = P f + AUa- 3 , 
P' = -duH 

= -2E t/t - 2U^V(x) - 2U^x d x V ^ 



<x 3 (P + <t 3 U^A) A - <r 3 [/t * dx (a x') 



where 



x' = Ucr 3 U n = Ua 3 P + U^UA (59) 

was used. In terms of coordinates and momenta, the constraint reads 

{Uct 3 P) 3 = . (60) 

Equation (|60|l is equivalent to © both in the presence and in the absence of a magnetic 
field. Taken together, l|58|l and (|59|l lead back to the equation of motion i|16[l . 
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Finally, let me mention an important subtlety regarding the component 
decomposition of the spinor equation. If, according to 0, U is represented as 
U = uq + Iu with u — Uk&k and pk denotes the momentum component 

conjugate to Uk, the spinor momentum is F = po — Ip with p = Y^k=iPk a k- 
The negative sign is necessary because in the spinor formulation the bivector IuPk 
is conjugate to IkUk- Dropping the bivector factors Ik leads to the stated result. 



4. The Kepler problem 

The unperturbed Kepler motion is described by the Hamiltonian 

H= -P^P - EU ] U = 2 , (61) 

The equations of motion derived from l|61|) are 

U' = pi, P' = 2EU f (62) 

or, as in JTSJ , 

U" = 2EU . 

If E < 0, this is the equation of motion of an isotropic four-dimensional harmonic 
oscillator, whose general solution reads 



U = A cos( \/-2E t)+B sin( -2E r) (63) 
with two constant even multivectors A and B bound, by 10 and (|61|) . to satisfy 

{Aa 3 B% = (64) 

and 

A^A + B^B = --. (65) 
E 

In the case of the pure Kepler motion, the angular momentum vector L and the 
Lenz vector e are conserved. Together, they uniquely specify an orbit HJ. I will now 
derive the KS-transformed expressions for these constants of motion. Throughout, the 
validity of the constraint will be assumed. 

The angular momentum vector is given by 

L = x x x = —I(xx)2 ■ (66) 

Within the geometric algebra, it is more convenient to introduce the angular 
momentum bivector 

1 = IL= (xx) 2 , (67) 

which specifies the orbital plane instead of the direction perpendicular to it. By 
© and ©, 

l = (^U* 3 tf U<t 3 U^) 2 
= \{UU% (68) 
= \{UP) 2 . 
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That / is conserved can be verified by a straightforward differentiation. 
Alternatively, it can be checked that the Poisson bracket vanishes, 

{/, H} = {dpH) * dul - (duH) * d P l 

= PUdul + 2EUUd P l 

1 ♦ + ( 69 ) 

= -(PiPh + E(U^Uh 

= . 

Note how the Poisson bracket formalism extends not only to multivector coordinates 
U and P, but also to non-scalar arguments. 
The Lenz vector is given by 



; • x 
e = Ix 

r 



(70) 



= IP^azlJ- 1 - Ua-sU- 1 . 
To calculate the Poisson bracket {e,Tt}, use {l,H.} — and 

PUd u U- 1 = -U- 1 P^U~ 1 (71) 

to find 

{e,H} = 1{P^(T 3 U- 1 ,H} - {U(T 3 U-\n} 

= l(PU d v (P^trsU- 1 ) + 2EUU dp (P^(7 3 U- 1 )) 

-PUdu (UctzU- 1 ) 
= ^ (-P t <T3C/- 1 pt C /-i + 2EU(t 3 U- 1 ) (72) 

- PVat/- 1 + E/oW^P+U -1 
= [1(-P^(t 3 U- 1 P^(t 3 U'' + 2EUU*) 

- Pt[/t + Ua 3 U- 1 P ] <T 3 U J{ jC/^Vaf/- 1 

Due to the constraint (|60|) . 

pt<7 3 [/t = Uct 3 P , (73) 
so that equation (|72|l simplifies to 

{e, H} = [l(-P*P + 2EU t U) - ptf/t + UP] U^a^U' 1 

= [-2lH + 2(UP) 2 ] U^cr-iU- 1 ( 74 ) 
= . 

Thus, the Lenz vector e is actually conserved. 

5. The Kustaanheimo-Stiefel description of closed orbits 

If the initial conditions £c(0) and x(0) for a trajectory are given, the pertinent initial 
conditions for the spinors U and V can usually be obtained, up to a choice of gauge, 
from © and (|13|l . This prescription fails for trajectories starting at the origin, where 
the KS transformation is singular. Although at first sight these trajectories appear 
rather exceptional, they bear a particular significance to atomic physics: Semiclassical 
closed-orbit theory [10111111X5] respresents an atomic photo-absorption spectrum as a 
sum over closed classical orbits, i.e. orbits that start at the nucleus and return to it. 
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Due to that particular importance, closed orbits deserve a detailed discussion. In this 
section, I will derive initial conditions for a trajectory starting at the Coulomb centre. 
I will then construct an orthonormal basis of the KS spinor space that is suitable 
for an investigation of the stability of orbits closed at the centre because it allows to 
separate the physically distinct directions along and transverse to the orbit from the 
gauge degree of freedom. 

To find initial conditions for orbits starting at the nucleus note that in the 
vicinity of the nucleus the Coulomb interaction is so strong that it dominates all 
external forces. The dynamics close to the nucleus is therefore described by the Kepler 
equation of motion l|18l) and its solution (|63(l . With the initial condition U(t = 0) = 
implemented, (|63|l reads 

U(t) = -^L= sm(V^2Er) . (75) 

Uq is the initial velocity in KS-coordinates. It must be normalized to 

U'M = 4 . (76) 

The choice of gauge for U' is arbitrary. 

The position vector corresponding to i|75|) is 

( \ 1 tt> tt<\ ^ 2 {V^Er) 
x(t) = ^^o^t/o — . (77) 

Thus, H75fl describes an electron moving out from the nucleus in the direction of the 
unit vector 

s = hufatf ■ (78) 

Uq is therefore a spinor rotating the vector <t% to the starting direction s and 
normalized according to (|76|l . In terms of the starting angles d and ip it is given 

by 

U' Q = 2 e -hv/2 e -h#/2 e -h«/2 (79) 

with an arbitrary gauge parameter a. The exponentials in l|79l) describe a sequence of 
three rotations taking the reference vector to the starting direction s. The initial 
momentum reads 

P =U$=2 e^ e h-9/2 e i 3V /2 (80) 

Its component decomposition is 



Pa = 


2 cos — 
2 


cos 


Pi = 


„ • 
2sm- 


sin ■ 


P2 = 


—2 sin — 
2 


cos 


P3 = 


—2 cos — 
2 


sin 





+ 


a 


t 


2 




<P 




a 




2 




Y 




a 




2 




Y 


+ 


a 



(81) 



2 

To describe the stability of a classical trajectory, a coordinate system with one 
coordinate along the trajectory and two coordinates perpendicular to it is customarily 
introduced in the neighbourhood of the trajectory. A linear stability analysis then 
requires calculating the derivatives of positions and momenta with respect to the 
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transverse initial conditions. Most conveniently, derivatives with respect to two 
orthonormal directions can be used. If these derivatives are to be calculated within 
the framework of the KS theory, for a given starting direction s and a direction s w _Ls, 
a KS spinor P u must be found such that a variation of the initial KS momentum Pq 
in the direction of P u corresponds to a variation of s in the direction of s w . 

As the initial momentum is given in terms of the starting angles in l)8()[l. the 
derivatives dPo /dd and 8Pq / dip can be expected to describe variations of the starting 
direction in the directions of increasing $ and ip, respectively. To check this and to 
find the correct normalization of the spinors, I will now construct, for a fixed direction 
s, a basis of the spinor space such that one of the basis spinors describes a variation 
of initial momentum along the orbit, i.e. in the direction of s, two give variations in 
the directions of two perpendicular unit vectors, and the fourth basis spinor describes 
the gauge degree of freedom introduced by the KS regularization. 

To this end, I consider a family of trajectories parameterized by an arbitrary 
parameter to. All trajectories start at the nucleus, the starting direction is given by a 
family of vectors s(uj). By (J7SJ), the initial KS momenta Pq(co) then satisfy 

a(o>) = -ptM^PoM , (82) 
so that in complete analogy with (jHJ) and Ijlljl 

0«_l/ pt dPc 
dto ~ 2 \ 3 dLU 



ds 1/ t dp \ 
sr. = o ( p o^rr7 ) ( 83 ) 



and 



dP Q 1 ds 1 ds 

a7 = r 3Po aJ ^ 2 PoS d^ (84) 



if the gauge condition 

P„Wg) s =0 (85, 

is imposed. Equation 1)84(1 gives the variation in the initial KS momentum pertinent 
to a given variation in the starting direction. If two different variations arc given, the 
scalar product of the momentum variations is 

«*>y. (86) 

duj\ J du>2 duji dt02 
so that the variations of KS momentum calculated from 184|) are orthonormal if the 
prescribed variations of the starting direction are. 
For a fixed starting direction 

s = e-^e-'^a^ 2 (87) 
given by the starting angles z9 and I now introduce the orthogonal vectors 

s v = e-^/V 2 e^/ 2 . (88) 

These are the unit vectors in the directions of ds / dd and ds / dip, respectively. The 
orthonormal basis s, s$ and s v of the position space gives rise to the three orthonormal 
KS spinors 

Pa =\czPoS = , 

=i<r 3 Po^= ^c-^Po , (89) 
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This set is complemented by a fourth orthonormal spinor 
F)P 1 

P a = ^ = -a 3 P I. (90) 

(Note the analogy with (JHSJ).) The spinor P a maximally violates the gauge condition 
(|85|l in the sense that 

PqWc = (P Wa>3 ■ (91) 

It therefore gives the direction in spinor space that corresponds to a gauge 
transformation, whereas P s describes a change of momentum along the orbit (i.e. 
a change of the energy) and P§ and P v give directions perpendicular to the trajectory. 
The desired separation of the physically distinct degrees of freedom has thus been 
achieved. Note that P$ = dP^/dd as anticipated, whereas dPo/dip = Po/3/2 does not 
satisfy Instead, 
dP 

= P Q costf + PoSini? . (92) 

dip 

In components, the four basis spinors read 
P s = cos I cos ^ - h sin § sin ^ + I 2 sin § cos ^ + h cos § sin 2+2. , 
P»=- sin I cos ^ - h cos | sin ^ + / 2 cos § cos ^ - 7 3 sin § sin ^ , 
P v = - sin I sin 2±£ - ^ cos f cos ^ - h cos f sin ^ + 7 3 sin f cos , 
P a = - cos f sin + Ja sin f cos J£=H + / 2 sin f sin ^ + h cos f cos *±2 . 

These formulae prescribe a basis of the spinor space uniquely up to the choice of a if 
?9 7^ 0, 7r. At the poles, the angle <p is undefined. Because in this case (|88[) gives a pair 
of orthonormal tangent vectors for any choice of ip, (|93|) can be used with arbitrary ip. 
The basis P s , Pg, P^, P Q of the KS momentum space can be supplemented by position 
spinors Uj — Pj to obtain the basis of a canonical coordinate system in spinor space. 
This basis set can then be used to investigate the stability properties of closed orbits. 
As a particular application, consider the stability determinant 

M' = det%^ (94) 

occurring in the closed-orbit theory description of atoms in crossed-fields |12I13I14I15| . 
It encodes the stability of a closed orbit starting from the nucleus in the direction 
characterized by the angles 1?^ , ip>i and returning from the direction As it stands, 

<|94[1 is not well suited to practical computations because it suffers from the singularities 
of the spherical coordinate chart: At the poles, neither the angle ip nor the angular 
momenta and p v are well defined, so that close to the poles, the calculation of the 
determinant becomes numerically unstable. With the help of the spinor basis H93|) . 
(|94|l can be rewritten in a form that is not plagued by any singularities. 

For a trajectory returning to the nucleus at time r — with a final KS momentum 
Pf, the solution l|75(l of the Kepler equation of motion takes the form 



p f 

U(t) = 



2E ■ ' V 2 f 



^2E y 1 V 2 ' ' (95) 

P( T ) = P f cos(V-2Pt) = Vl + ErP f , 

which is valid as soon as the electron is sufficiently close to the nucleus so that the 
external fields can be neglected. In l|95|l . the normalization condition U^U — 2r and 
the pseudo-energy conservation l|61(l were used. 
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To discuss the stability properties of an orbit, consider a family of orbits 
parameterized by an arbitrary parameter u>i . Equations of motion for the derivatives 
dU /Olui and dP/ duj\ are then obtained by linearizing the equations of motion for 
U and P around the unperturbed trajectory. Close to the nucleus the dynamics in 
governed by the linear equations l|62[l . so that the perturbations satisfy the same 
equations of motion as the coordinates and momenta. In analogy with ()63[) . they read 

dU _ M\ 



duJi ^2E 



sin( v / -2£' t) + M 2 cos(y/-2E r) , 



dP 

— = Mi cos{^/-2E t) - V-2E m\ sin(V ' -2E) 

OU}\ 



(96) 



with constant even multivectors Mi and M2. With the help of the normalizations 
inferred from 1)95(1 . equation (|96Jl simplifies to 

— =-,/- Ml + VT+EiM 2 , 

dw-L V 2 ^ 
= ^1 + Er Mi - V2rE M\ 

OLUl 

so that the constants Mi = dPf/duji and Mi = dUf/du)\ can finally be identified 
with the values of the derivatives obtained at r = 0. 

According to l|68|) . the angular momentum component in a plane specified by a 
bivector B is 

L B =L*B = ^(BUP) . (98) 

The derivative of Lb with respect to a parameter u>i can then be calculated with the 
help of (|9"7|) . It is given by 

It does not depend on r because all components of the angular momentum bivector 
are conserved in the Coulomb region. 

For the angular momentum component p lxJ2 conjugate to an angular coordinate 
0J2, the relevant bivector is 

B = sp-, (100) 



duo 



2 



a v „ i/ P *«!i\_/ P n h (101) 



so that 

dui 2 \ 1 duj 2 dcui I \ W2 dui 

where P u , 2 denotes the basis spinor corresponding to 0J2 by H84|) with the final 
momentum Pf used in place of Pq . 

Consider coordinates ■& and <p, such that ds/dd = 84 and ds/d(p = s v are unit 
vectors given by l|88|l . The stability determinant M' can then be rewritten as 



sin di 



det " ? T = det % » = smtiiSintifM (102) 
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with a 2 x 2-determinant 




free of any coordinate-induced singularities. It is this form of the stability determinant 
that was used in |T5] . 

The derivatives dUf/du) needed in i|103|) can be calculated numerically by 
integrating the linearized equations of motion 

d du (du \ (dp \ 

-j — tt - = [-^- * du) dpH + — * dp) dpH , 
dr ouj \ouj / \ouj I 

v 7 v 7 (104) 

d dP fdU \ fdP \ 1 j 

dr V aw / V ouj I 



along the closed orbit with initial conditions 



(0)=0, ^_( ) = P W . (105) 



6. Conclusion 



The geometric-algebra formulation of the KS transformation endows the four KS 
coordinates with a clear geometric interpretation. It was shown in this paper that it 
also allows for an easy incorporation into a Lagrangian and Hamiltonian formulation 
of the KS theory. At the same time, the known KS Hamiltonian, that was restricted to 
homogeneous external fields, was generalized to describe the KS motion in arbitrary 
static external electromagnetic fields. This result is of interest beyond the realm 
of atomic physics because it is equally applicable to an arbitrary conservative non- 
electromagnetic force. 

Closed orbits starting at and returning to the Coulomb centre require a special 
treatment that takes the singularity of the KS coordinate frame into account. The 
geometric algebra offers a convenient way to derive closed orbit initial conditions and 
a basis of the spinor space that separates physically distinct degrees of freedom. 

The calculations carried out here amply illustrate the power of the geometric 
algebra formulation of the KS theory. It can thus be expected to form a convenient 
starting point for analytic calculations in classical perturbation theory. 



Appendix A. Introduction to geometric algebra 

Geometric algebra is an algebraic system designed to represent the geometric 
properties of Euclidean space in the most comprehensive and systematic way possible. 
It was pioneered by Hermann Grassmann and William Kingdon Clifford during the 
nineteenth century. From the 1960's on, David Hestenes, with the aim of providing 
a universal mathematical framework for theoretical physics, extended the algebraic 
techniques of Grassmann and Clifford by a differential and integral calculus within 
the geometric algebra, which he called geometric calculus |16j . 

This appendix gives only a sketch of geometric algebra as far as it is needed in 
the present work. A more extensive introduction, with an extension to Minkowski 
spacetime, is contained in 18, 17 . A thorough introduction to the geometric algebra 
of Euclidean 3-space, with a detailed discussion of applications to classical mechanics, 



The Kustaanheimo-Stiefel transformation in geometric algebra 



18 



can be found in 0. A detailed presentation of the mathematical properties of the 
geometric algebra is given in |16j . 

Appendix A.l. The geometric algebra of Euclidean 3-space 

The orientation of two vectors a and b in space can be characterized by the projection 
of one vector onto the other, which is described by the scalar product a ■ b, and the 
plane spanned by a and b, which is characterized by the vector product a x b. In 
the geometric algebra these complementary products a ■ b and a x b are unified into a 
single "geometric" product ab. The construction of the geometric product starts by 
picking a right-handed frame of orthonormal unit vectors eri, <t 2 and cr 3 . For them, 
the existence of an associative, but non-commutative geometric product satisfying 

rr rr . + rr.cr, = 25 i3 (A.l) 

is assumed. In addition, the geometric product is required to obey the distributive 
law with respect to the usual addition of vectors. It follows from 1A.1|) that <xf = 1 is 
a scalar. The reader may notice that the defining relation (|A. If) is the same as obeyed 
by the Pauli spin matrices. Indeed, these matrices generate a matrix representation 
of the Clifford algebra of Euclidean 3-space. In the present context, however, it 
is important to retain the interpretation of the er.j as ordinary vectors instead of 
regarding them as matrices. The elements of the Clifford algebra are thus given a 
geometric interpretation, as is indicated by the name "geometric algebra" introduced 
by Clifford himself. It turns out that all calculations within geometric algebra can be 
done without recourse to a matrix representation. 

By virtue of the defining relation (|A.1|I . the geometric product of two arbitrary 
vectors a = J~li—i GjtTj and b — y^._ 1 bi<Ji is 

ab = ai&i + a 2 b 2 + a 3 b 3 (A. 2) 

(a 2 6 3 - a 3 b2)a-2(T3 + (a 3 bi - aib 3 )cr 3 (Ti + (aib 2 - a 2 &i)<T 1 er 2 . 

The scalar terms of this equation comprise the scalar product a ■ b. In addition, 
there are terms containing the product of two orthogonal vectors. These terms are 
neither scalars nor vectors. They are referred to as bivectors. As their coefficients 
are the components of the vector cross product a x b, bivectors should be interpreted 
as describing an oriented area in the same way as a vector describes an oriented line 
segment. Accordingly, a product of three orthonormal vectors is called a trivector and 
interpreted as representing an oriented volume element. With the help of the unit 
trivector 

/ = erier 2 <T 3 , (A.3) 

equation (|A.2fl can be rewritten as 

ab = a ■ b + la x b , (A.4) 

which achieves the desired unification of the scalar and vector products. Notice 
that (|A.4Jl contains a sum of quantities of different types, a scalar and a bivector. 
This should be regarded as a formal sum combining quantities of different types into 
a single object with a scalar and a bivector part, in analogy to how a real and an 
imaginary number are added to yield a complex number. 
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As the scalar product is symmetric in its factors whereas the vector product is 
anti-symmetric, these products can be recovered from the geometric product via 

a b = — (ab + ba) , 

\ (A.5) 
a x b = — (ab — ba) . 

In particular, parallel vectors commute under the geometric product whereas 
perpendicular vectors anti-commute, and any vector a satisfies aa = a ■ a. 

It is a crucial feature of the geometric algebra that it contains elements of different 
grades, viz. scalars, vectors, bivectors, and trivectors. A general element can be 
written as a sum of these pure-grade components and is called a multivector. The 
pure-grade parts of any multivector can be extracted by means of a grade projector. 
Let (A)k be the grade-fc part of the multivector A. Due to its particular importance, 
the scalar projector can be abbreviated as (A) — (A)o, and the scalar product of two 
multivectors A and B is defined by 

A*B = (AB) . (A.6) 

For vectors, this agrees with the scalar dot product. Any two multivectors commute 
under the scalar product: 

A*B = B* A . (A.7) 

In a term containing different kinds of products, the scalar product as well as the 
vector cross product are understood to take precedence over the geometric product. 
This convention has already been used in (|A.4|I . 

A multivector which contains only parts of even grades, i.e., scalars and bivectors, 
is referred to as an even multivector. The even multivectors form a subalgebra of the 
full geometric algebra. 

Finally, the reversion A^ of a multivector A is obtained by interchanging the order 
of vectors in any geometric product. Thus, bivectors and trivectors change sign under 
reversion, whereas scalars and vectors remain unchanged. Formally, the reversion can 
be defined by the properties a^ = a for any vector a and 

(AB)* = B*A* , (A + B^=A^+B^ (A.8) 

for multivectors A and B. 

Appendix A. 2. Rotations in the geometric algebra 

Within the geometric algebra rotations are conveniently represented in the form 

a i ► R(a) = RaR) (A.9) 

with an even multivector R satisfying the normalization condition 

= 1 . (A.10) 

Conversely, any normalized even multivector describes a rotation. An arbitrary even 
multivector satisfies a — UU^ > 0, so that U = yfaR is a multiple of a rotor R. 
Therefore, 

UaU^ = aRaR^ , (A.ll) 
and U describes a rotation-dilatation of 3-space. In particular, 

UaU* = (Uarf)i (A.12) 
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is a vector for any even multivector U and any vector a. 

A rotation can be characterized by specifying two vectors a and b so that a is 
mapped to b by a rotation in the plane (ab) 2 spanned by a and b. The rotor R 
describing this rotation is 

R=^= / 1 + ba ■ (A.13) 
\a + b\ y/2(l + a-b) 

If, alternatively, a rotation is characterized by its rotation axis, given by a unit vector 
n, and a rotation angle ip, the pertinent rotor reads 

R = c -i nv /2^ ( A14 ) 

where the exponential function of an arbitrary multivector is defined by the familiar 
power series 

v An 

c- 4 = £^- (A.15) 
It satisfies the "power law" relation 

e A+B = e A e B (A 16) 

if AB = BA and 

e A B = Be A if AB = BA , (A.17) 

e A B = Be~ A if AB = -BA . (A. 18) 

Appendix A. 3. The multivector derivative 

The formalism of the multivector derivative provides a differential calculus for 
arbitrary multivector functions. Let F(X) be a smooth multivector- valued function of 
the multivector argument X. Neither the grades contained in X nor in F are specified. 
The directional derivative in the direction of a fixed multivector A by 

dF(X + rP x (A)) 



A * d x F(X) 



dr 



(A.19) 



where Px(A) projects A onto the grades contained in X. 1A.19|) agrees with the 
familiar definition of the directional derivative. 

Let ej, J = 1, . . . , 8 be a basis of the geometric algebra and e J its dual basis, i.e., 
ej * e K = 5*f . The multivector derivative is then defined to be 

d x = Y,e J ej*d x . (A.20) 
j 

It inherits the algebraic properties of its argument X. In particular, dx contains the 
same grades as X. Notice that the scalar product A * dx is indeed the directional 
derivative in the direction A, justifying the notation introduced in (|A.19() . For a vector 
argument x, the multivector derivative d x reduces to the vector derivative, which is 
analogous to the familiar nabla operator. 

Both the directional derivative and the multivector derivative are linear operators 
and satisfy Leibniz' rule 

A*d x (F(X)G(X)) = (A* d x F(X))G{X) + F(X)(A* d x G(X)) , (A.21) 
d x (F(X)G(X)) =dxF (X)G(X)+ *dx F(X) G (X) . (A.22) 
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In HA.22[) . the overstars indicate the functions to be differentiated. Notice that the 
second term in (|A.22|I is in general different from F(X)(dxG(X)^j , because due 
to its multivector properties the multivector derivative does not commute with F 
even if F(X) is not differentiated. The directional derivative, on the contrary, is 
a scalar differential operator that commutes with any multivector that is not to be 
differentiated. For this reason it is often convenient to write the multivector derivative 
as 

d x = d A A * d x ■ (A.23) 

This form decomposes dx into a multivector &a and a scalar differential operator 
A * dx, which can be moved freely among multivectors. 

In addition, the directional derivative satisfies the chain rule 

A * d x F(G(X)) = (A * d x G{X)) * d G F{G) , (A.24) 

which is useful in many calculations. 

A fundamental result concerning the multivector derivative is 

d x (XA) = d x (AX) = P x (A) (A.25) 

for any multivector A. As a consequence, 

dx(X^A) = dx(AX^) = Px(A^) . (A.26) 
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